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Abstract 

Consistently with the Einstein equivalence principle and using an electromagnetic 
formulation first suggested by Tamm, we show that a local observer cannot observe 
negative-phase-velocity (NPV) propagation of electromagnetic waves in vacuum, whereas 
a global observer can appreciate that phenomenon. Using the specific example of the 
Kerr metric, we also demonstrate the possibility of NPV propagation within the ergo- 
sphere of a rotating black hole. 
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1 Introduction 

Recently, we have shown that electromagnetic plane waves can propagate in classical vacuum 
in such a way that the phase velocity vector has a negative projection on the time-averaged 
Poynting vector, provided that the vacuum is nontrivially affected by a gravitational field [1]. 
Our approach is based upon Tamm's electromagnetic formulation [2], involving a gravito- 
magnetic metric. The negative-phase-velocity (NPV) propagation characteristic depends 
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on the direction of the propagation wavevector but not on the frequency. NPV propagation 
characteristics in certain homogeneous material mediums lead to the phenomenon of negative 
refraction [3]-[5], which suggests the possibility of gravitationally assisted negative refraction 
of electromagnetic waves by vacuum, with attendant scientific implications [1],[6]. 

Researchers have studied electromagnetic wave propagation in terrestrial environments 
for several centuries, and our planet is certainly affected by the solar gravitational field. Yet, 
NPV propagation in classical vacuum has never been reported in the terrestrial context. 
Although the solar gravitational field evidently does not satisfy the conditions for NPV 
propagation that were reported elsewhere [1], evidence accumulated from the Pioneer 10 
mission has prompted suggestions that new physics may be found even at the outer reaches 
of our solar system [7] . These and other considerations beg the question: Are global and local 
perspectives of gravitationally assisted NPV propagation different? In this communication, 
we answer that question in the affirmative. 



2 Maxwell equations in universal spacetime 

In the absence of charges and currents, electromagnetic fields obey the covariant Maxwell 
equations^ 

/q/3;i/ + //3i/;a + = , /i"^.^ = 0, (l) 

where fa/s and are, respectively, the covariant and the contravariant electromagnetic field 

tensors whereas the subscript indicates the covariant derivative with respect to the z/th 
spacetime coordinate. The spacetime — denoted by the vector containing the normalized 
time coordinate = ct (with c as the maximum speed of light in the absence of the 
gravitational field) and the space coordinates x^, x^, and x^ — is Riemannian, with the 
metric gai3 being a function of spacetime and carrying the signature (+,—,—,—) [8]. 

It is commonplace to follow up on a suggestion of Tamm [2],[9]-[ll] and change the form 
of (1) for application to electromagnetic fields in vacuum. The Maxwell equations (1) may 
be expressed in noncovariant form as 

fap,iy + fpu,a + fva,l3 — , 

wherein g — det [gai3\ and the subscript denotes ordinary differentiation with respect 
to the vih spacetime coordinate. We note that generalizing the Maxwell equations from 
noncovariant to covariant formulations is not totally unambiguous [12], [13]. In the absence 
of experimental evidence to eliminate this ambiguity, we adopt the standard generalization 

(!)• 

Let us introduce the electromagnetic field vectors E"^, S^, and i?^ via the standard 
decompositions 

= fio , Bj^ — {l/2)eirnnfmn 

Dj^i-gf'h'' , H, = (l/2)£w i-gf' U 

§ Roman indexes take the values 1, 2 and 3; Greek indexes take the values 0, 1, 2, and 3; summation is 
implied over any repeated index; and Gaussian units are used. 



= 0, 



(2) 
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with Ejimn being the three-dimensional Levi-Civita symboL Thereby the noncovariant Maxwell 
equations (2) assume the familiar form 



Bt,l — , B(,fi + eirnnEm,n — 
Di^e. — , — -D^,o + ^emnHm,n — 



(4) 



In vacuum, the components of the electromagnetic field tensors are connected by the 
constitutive relations 

Uf.^g^.gf.^h'^''. (5) 
These constitutive relations of vacuum can be stated for the electromagnetic field vectors as 



Hn 

Be = 'yimHm — ^emn Tm En 

where 



(6) 
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iim = - {-gf^ — , (7) 

S'oo 

= (8) 

Equations (1) and (5) employ curved spacetime. So do (4) and (6), but the difference 
is that they look like the familiar electromagnetic equations in flat spacetime applied to an 
instantaneously reacting, bianisotropic medium. Techniques commonly employed to handle 
electromagnetic problems in the absence of gravitational fields should therefore be useful for 
solving (4) and (6). Before proceeding further with our analysis, let us therefore recast (4) 
and (6) using the conventional 3- vectors and 3x3 dyadics as 



and 



d 

V X E(ct,r) + TrBict,r) = 

at 

d 

V X H{ct,r) - -^D{ct,r) = 



D(ct,r) = eoj(ct,r) ■ E(ct,r) - -T(ct,r) x H(ct,r) 
= c 

B(ct, r) = iio l(ct, r) ■ H(ct, r) + - r(ct, r) x E(ct, r) 
= c 



(9) 



(10) 



wherein j{ct,r) is the dyadic-equivalent of 7^^, T{ct,r) is the vector-equivalent of and 
SI units are adopted. 

3 Global perspective 

Suppose that we wish to solve these equations in a certain region X of spacetime, subject 
to specific boundary conditions. A fairly standard procedure would be to partition X into 
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subregions ("^^A", {n — 1, 2, 3, ... ), in each of which we would replace the nonuniform metric 
Qa^ by the uniform metric ^""^gaB- Correspondingly, the nonuniform vector with components 

would be replaced by the uniform vector with components ("^F^ in each subregion ("^A". 
The spacetime a;" would still remain curved. After solving (4) and (6) in each subregion, 
we could stitch back the subregional solutions into the regional solution. This piecewise 
uniform approximation technique is very common for solving differential equations with 
nonhomogeneous coefficients [14]. 

Accordingly, let us examine the electromagnetic fields in the nth subregion. A three- 
dimensional Fourier transform of the electromagnetic field vectors can be taken, with the 
wavevector k denoting the Fourier variable corresponding to r; thus, 

-j^ POO poo poo 

E(a;/c, k) exp [i{k"r — cut)] duj dki dk2 



E{ct,r) = 



POO POO POO 

H_{ct, r) — - / / H(<^/c, k) exp [i{k'r— ujt)] duj dki dk2 

^ J —oo J —oo J —oo , 



(11) 



where % = and cu is the usual temporal frequency. The wavevector component k^ is 

determined by first substituting (11) in (9), and then formulating a 4x4 matrix ordinary 
differential equation which is then solved as an eigenvalue problem [15]. 

As our focus lies here on propagating plane waves, evanescent (nonuniform) solutions 
are eliminated from further consideration by enforcing the restriction k^^ G M. We remark 
that E(a; / c, k) and H{u! / c, k) are the complex- valued phasors of the electric and magnetic 
fields, respectively, associated with a plane wave, for which a time-averaged Poynting vector 
{P_{uj/c,k))t can be derived. When the projection of k on {P.{uj/c,k))t is negative, i.e., 
k • {P_{oj/c^k) )t < 0, we say that the phase velocity is negative. 

The crucial quantity k_ • {P,{uj/c, k))t is derived as follows [1]: Combining the Fourier 
representations (11) with the constitutive relations (10) and the Maxwell curl postulates (9), 
we get 

(-)7.H(c./c,A:), 



p X E{Lj/c,k) 
P X H(a;/c,^) 



uJHo 



where 



p — K 



k--^-^r. 



(12) 

(13) 
(14) 

(15) 



The use of (12) to eliminate H(a;/c, ^) from (13) provides, after some manipulation, the 
eigenvector equation 



V c 



P 



1-P 



l + pp .H^l . E(o;/c,^) = 



for representing E(a;/c, ^), and the corresponding dispersion relation 



p • ^"^7 'p - ^— j 



0. 



(16) 



(17) 
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Herein, 7 is the identity dyadic, and |*^")7| is the determinant of ^"^7. Considering (16) in 
hght of (17), we see that all E(a;/c, ^) eigenvectors solutions must satisfy the relation 



P 



.(")7 • E(a;/c,fc) = 0. 



(18) 



Let us introduce the eigenvalues ^"'^71,2,3 and corresponding eigenvectors ^"^^1^2,3 ^"^7- 
For the purposes of illustration, suppose we consider planewave propagation along the direc- 
tion parallel to *^"^^3, while the vector ("^F lies at an angle 6 in the plane of the eigenvectors 

^'^'>v_i and '•"'^53, as per 

(n)f^ = (»)f sin^ ^ 

(")f2 = > . (19) 

Hfa = Hf cos^ ^ 

After exploiting the orthogonality condition (18), the general solution to (16) may be ex- 
pressed as the sum of two independent modes as 



E(a;/c,^) = Aa{uj/c,k)e^{(jj/c,k) + Af,{uj / c,k) ej^{uj / c,k) 



where 



e,{uj/c,k)= 



(1 1^1 - ^ Hf cos^) + ("^71 ^ (")f sin^ ("^Cg 



(20) 



(21) 



The complex-valued scalars Aa^ are unknown amplitude functions that can be determined 
from initial and boundary conditions. The corresponding general solution for H.(^/c, ^) 
follows straightforwardly by combining (20) with the Maxwell curl postulates (9) as 



H.(a;/c,^) — Aa{u; / c, k) h^{u; / c, k) + Ah{u; / c, k) hf,{u; / c, k) 



(22) 



with 

ha{^/c,k) = 

hii^ /c,k) = 



1 



UJHo 

1 



(^^ (")f cos^ - 1^1) - ^ (")f sin^ ("^^3] 

_ ^ Hfcos^)' ('^)73+ (")fsin^)' (" 



^7i 



(23) 

The dispersion equation (17) reduces to a |fc|-quadratic expression which yields the two 
wavenumbers 



|A;| 



UJ 



( 



(n) 



(")f cos^ ± W (")7i H72 - HPsin^e 



V 



(n) 



73 



(24) 
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Finally, combining (20), (22) with (24), we find that 

1 



k'{Piulc,k))t = 



(n)~, (n)~ 
7l 72 



(n)^3 



± ("^fcos^W W71 (") 



72 



(")r2sin2^ 



(n) 



73 



X (|^(a;/c, fc)|2 + ^736.^ j^jl \A,{uj/cM^) . (25) 

Let us probe the physical significance of (25). To do so, we focus upon a specific gravito- 
magnetic metric, namely the Kerr metric. We begin by noting that for a rotating black hole 
of geometric mass rurbh and angular velocity arhh in conventional units [16], the Kerr metric 
yields 



(26) 



and 



("^71= ^"^72 
W73 = 1 



+ 



(27) 



where '^"'^a;, '^"^y and ^'''•'^z are the values of the three space coordinates at some representative 
point in the subregion ("•'A'. The quantity '^"•'i? is defined implicitly via 



{n)^2_ (n)~2^ (n) ~2 ^ (n)~2 



(28) 



Since both |*^")7 | > and '■"'-'73 > 0, we see from (25) that NPV propagation in ("•'A' arises 
provided that the inequality 



I _ Hf . W;y . (")f <; 
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(29) 



is satisfied. 

It is illuminating to recast (29) in terms of the metric components as per 



(n)~00 



(n) 



goo 



< 0. 



(30) 



For the Kerr metric we have [16] 

(n)~00 ^ (n)^4 ^ ^2^^ {n) ~2 ^ 2mrbh ^"^^^ 
(")^00 ~ + a2(,^ _ 2^^^^ (n)^3- 



(31) 
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Within the ergosphere — which hes between the outer event horizon and the stationary hmit 
surface — we have 



and > a^^^^. Therefore, inside the ergosphere 

(")^ > 

and the NPV condition (30) is evidently satisfied for the chosen wave-propagation case. 

Thus, when the traversal of an electromagnetic signal is traced from point P E ^^^X to 
point Q e ^'^^X, with p q, the possibility of NPV propagation in some subregions of X 
cannot be ruled out a priori. 




4 Local perspective 

The previous two sections are from a global perspective. Let us now attend to the local per- 
spective of an observer located at some point p. This observer is constrained to formulate 
an electromagnetic theory valid only in some small neighborhood of p, and this theory must 
emerge from observations made only in that neighborhood. The gravitational field in this 
neighborhood can be held to be virtually uniform, the metric from the global perspective 
being gadlp- However, given the admissibility of the uniform-gravity theory for all local ob- 
servers in this neighborhood, we show here that the local observer would end up formulating 
a local metric — which is the same as in the special theory of relativity. 

By virtue of the Einstein equivalence principle [17], ga/slp is constrained to be such that 
there exists a matrix yielding 

?7^. = A^^a/3U (34) 

where rj^^, = diag [1, —1, —1, —1] is the Lorentzian spacetime metric [8]. The matrix A^ can 
be used to construct a flat spacetime x'" for all x" in the specific neighborhood of p as per 



^ ~ dxP 

dx^dx^dx^dx^ = det [A^] dx'^dx'^dx'^dx''"^ 



(35) 



We must note that, as A" depends upon ga/slp and therefore on X I p , ci global coordinate 
transformation which simultaneously transforms Qa/s into rj^iy for all x" E X cannot be 
realized unless gravity is totally ignored [8]. 
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In the local coordinate system x'", the spacetime metric in the specific neighborhood is 

given as rja/s- Repeating the steps outlined in Section 2 with r]ai3 substituting for g^js, we find 
that the local constitutive relations for the neighborhood reduce to the trivial form 



with respect to x'°'. NPV propagation is not possible for the medium characterized by the 
constitutive relations (36) [18]. 

5 Concluding remarks 

Thus, we have shown that gravitationally assisted negative-phase-velocity propagation in 
vacuum can be appreciated only from a global perspective based on curved spacetime; 
whereas a local observer, constrained to a flat spacetime, must conclude NPV propaga- 
tion is impossible in vacuum. The diversity of the local and the global perspectives is in 
full accord with the Einstein equivalence principle, and can be explained as follows: The 
subregional metric ^'^^Qap appears from a suitable subregional averaging of Qa^. Whereas 
must satisfy the dictates of the Einstein equivalence principle at every e X, there is no 
reason for ^'""^^Qap to satisfy the same dictates at even one x°' E ^"^^X. 

The feasibility of NPV propagation within the ergosphere of a rotating black hole has 
been demonstrated in this communication. A detailed numerical study is called for in order 
to further explore the circumstances and implications of such type of propagation, which we 
hope to report in due course of time. 
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